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WILSON, J. W., and TOWNSEND, L. W. A Benchmark for Galactic Cosmic Ray 
Transport Codes. Radiat. Res. 

A nontrivial analytic benchmark solution for galactic cosmic ray 
transport is presented for use in transport code validation. Computational 
accuracy for a previously-developed cosmic ray transport code is established 
to within one percent by comparison with this exact benchmark. Hence, 
solution accuracy for the transport problem is mainly limited by inaccuracies 
in the input spectra, input interaction databases, and the use of a straight 
ahead/vel oci ty-conservi ng approxi mati on . 
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INTRODUCTION 

With the advent of the Space Station era, and future manned lunar bases 
and Mars missions under consideration, considerable attention must be given to 
developing methods for shielding against the high-energy heavy ion (HZE) 
component of galactic cosmic rays, since their ranges are generally comparable 
to their mean free paths for nuclear interaction in tissue and low density 
shield materials. Because these HZE particles include nearly all nuclear 
species, and possess a broad spectrum of energies, detailed laboratory 
measurements for all possible ion-shielding combinations are not practical; 
therefore, accurate calculational methods describing the interactions and 
transport of these energetic ions in bulk matter are needed to properly 
evaluate the shielding effectiveness of spacecraft structures and the self- 
shielding factors of the astronauts themselves. 

Whenever galactic cosmic rays traverse bulk matter, their radiation 
fields change composition through interactions with the target materials 
encountered. Aside from continuously losing energy through collisions with 
orbital electrons, the incident ions and target nuclei also undergo nuclear 
fragmentation (breakup) reactions. These fragmentations result in the 
production of secondary and subsequent generation reaction products which 
alter the isotopic composition of the transported radiation field. Studies of 
these radiation fields are presently hampered by the lack of an adequate 
nuclear fragmentation cross section data base over the broad spectrum of 
energies and fragmenting species involved. Because the experimental 
fragmentation data base is sparse and no adequate quantitative fundamental 
theory exists, a semiempirical formulation was developed for use in cosmic ray 
studies (1). Unfortunately, it lacks charge and mass conservation for 
secondary fragments (2). Recently, an alternative semiempirical fragmentation 
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model, which is based upon more fundamental physical considerations and which 
conserves fragment mass and charge, has been formulated (3). Its predicted 
cross sections agree with available experimental fragmentation data to the 
extent that these data agree among themselves. This fragmentation model has 
been successfully incorporated into a recently-developed galactic cosmic ray 
(GCR) transport code (4,5) for use in space radiation shielding 
applications. In the GCR code, methods previously developed for nucleon 
transport (6) were extended to HZE transport by a combination of analytic and 
numerical tools. The GCR ion transport problem was transformed to an integral 
along the characteristic curve of that particular ion. Fragment velocity 
conservation then enables the perturbation series (6) to be replaced by a 
simple numerical procedure. The resulting method reduces the difficulty 
associated with low-energy discretization, eliminates any restrictions to a 
particular functional form for the stopping power, and is computationally 
simple and nondemanding upon computer resources. Details of the solution 
method and a comparison with an alternative code are published elsewhere 
(4,5). 

In the present report we address the question of GCR transport code 
validation. Ideally, validation should be accomplished using detailed 
transport data obtained from carefully planned and controlled experiments; 
unfortunately, there exists a paucity of such data. Although useful for 
comparison purposes, the atmospheric propagation measurements used previously 
(5) are clearly not definitive since they consist of integral fluences of as 
many as ten different nuclear species combined into a single datum. Although 
limited quantities of HZE dosimetry measurements from manned space missions 
(e.g. Skylab) are also available (7), numerous assumptions concerning the 
relationships between dosimeter locations and spacecraft shield thicknesses 
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and geometry must be made in order to estimate astronaut doses using 6CR 
codes. Since many of these assumptions may involve inherently large 
uncertainties (a factor of two or greater), it becomes difficult to attribute 
sources of any comparison differences to particular assumptions or 
approximations which may have been used in the analyses. In the absence of 
definitive GCR transport measurements with which to compare code predictions, 
other methods of validation must be considered. As noted in references (2) 
and (4), there are several different versions of HZE transport codes 
available. When used with the same input spectra, interaction parameters, 
and boundary conditions, all should yield comparable results. The history of 
transport code development, however, suggest otherwise. For this reason, a 
realistic, nontrivial, exact, analytic solution to the simplified Boltzmann 
equation used to describe HZE transport has been formulated as an absolute 
standard for code comparison purposes. 
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GCR TRANSPORT BENCHMARK 


In passing through bulk matter, heavy ions lose energy through 
interactions with atomic orbital electrons along their trajectories. On 
occasion there is a violent collision with nuclei of the target medium. These 
collisions produce projectile fragments moving in the forward direction and 
low-energy fragments of the struck target nucleus which are nearly 
isotropically distributed. The transport equations for these short-range 
target fragments can be solved in closed form in terms of collision density 
(8); therefore, projectile fragment transport is the main subject of current 
interest in HZE transport. 

For the GCR transport problem, one typically uses the straight ahead 
approximation and neglects target secondary fragments enabling the transport 
equation to be written as (4,5) 





(i) 


where $.(x,E) is the flux of ions of type j with atomic mass A* at x moving 
J J 

along the x axis at energy E in units of MeV/amu, c. is the corresponding 

J 

macroscopic nuclear absorption cross section, Sj ( E ) is the change in E per 
unit distance, and m^ is the multiplicity of ion j produced in collision by 
ion k. The range of the ion is given as 


R. ( E) = / 
J 0 


dE 1 

S.(E') 


( 2 ) 



The solution to Eq. (1) is found to be subject to boundary specification at 
x = 0 and arbitrary E as 


$j(0,E) = Fj(E) 


where Fj(E) is called the incident beam spectrum. 
It follows from Bethe's theory that 


AZ* - 

S.(E) = -N— S (E) 

3 A.Z 2 p 
JP 


(3) 


(4) 


which holds for all energies above 10 MeV/amu provided the ions remain fully 
stripped. We can then write for the ji!l ion 



(5) 


where the subscript p refers to protons. Equation (5) is accurate for high 
energies but only approximately correct at low energy because of : (a) 
electron capture by the ion which effectively reduces its charge, (b) higher 
order Born corrections to Bethe's theory, and (c) nuclear stopping at the 
lowest energies. The range scale parameters v. are obtained from 

J 
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and are generally energy-dependent. When the ion velocity is large compared 
to the velocity of the orbital electrons the v. approach 

J 


For the benchmark problem, the incident spectrum is limited to a single 
ion type (j = J). Since the GCR spectrum for a typical ion is of the form 

F(E) ~ E" a (8) 


where a «» 2.5, we choose the energy spectrum to be of similar functional form 
as 


Fj(E) = «jj/(t R j(E)3 2 S j ( E ) } (9) 

Defining the characteristic variables 

nj = x - Rj(E) (10) 


and 


= x + Rj(E) (11) 

equation (1) can be solved by the method of characteristics (4,5) to give 
Tj(x,E) = e ajX /[vjX + Rj(E)] 2 


(12) 



where 


4>j ( x , E ) = Sj (E) $j(x,E) 


(13) 


and 


a J " 


(14) 


This is the trivial solution for the incident beam species. For j<J 
(secondary fragments) it can be shown that 


A A 

_ v i "( 0 i r 'i + SO/2 

*,(x,E) = o-iM • , iJx.E) e 3 J d 3 

J J JJ Vj J 


(15) 


where, in terms of the exponential integral function E 2 (x), 


I j( x,E) 


e*b lv j* VjUj/Z E 2 (bv ltl ) E ? (bv jSj ) 

( V J‘ v j> v j 5 j Vj 


(16) 


for j = J-l and 


b = 


,v j- V 


(17) 


Clearly, equations (16) and (17) are true for all j if m^j = 0 for all j<J 
(i.e. the secondary fragments themselves do not fragment). 
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BENCHMARK RESULTS 


The benchmark solution was calculated for an incident iron beam (J = 26) 
in an aluminum target, for which the input parameters are 


A ■- 

a 26 = cm/g 

2 

°25 = * 84 260 cm/g 

2 

M 25,26 °26 = - 00403 cm ^ 9 


Results of the GCR transport code simulation of this benchmark for the 
propagating incident iron beam and secondary manganese (j = 25) ions are 
displayed in Tables 1 and 2 where they are compared to the exact analytic 
predictions obtained from equations (12) and (15). It is clear from these 
tabulated results that the numerical solution methods developed previously 
(4,5) are accurate in solving equation (1) for GCR transport to within one 
percent. This indicates that any limitations to accurately solving GCR 
transport problems must focus upon the simplifying approximations used to 
obtain equation (1), as well as upon unresolved issues concerning the need to 
include multi pi e-Coulomb scattering effects, fragment momentum dispersion 
effects, and perhaps most importantly, the nature and quality of the input 
cross section data bases. To illustrate this last point, we are aware of only 
one heavy ion transport code (2) which uses energy-dependent cross sections. 
Recent studies, however, suggest that fully energy-dependent cross sections 
may be important for some transport code applications (9). 
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CONCLUDING REMARKS 


The need to develop suitable benchmarks for use in validating and 
comparing existing galactic cosmic ray transport codes has been described and 
an exact nontrivial analytic benchmark solution presented. This benchmark 
solution was then used to establish computational accuracy for a previously- 
published cosmic ray transport code to within one percent. Finally, remaining 
unresolved issues in GCR transport were briefly described. 
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TABLE 1 

Comparison of the Benchmark Numerical Simulation to the Analytic Solution for 
Iron Ions as a Function of Ion Energy and Depth into the Aluminum Abosrber. 



* Fe (0 ’ E) 

*p e (10g/cm?E) 

$ Fe (20g/cm?E) 

E,MeV/amu 

numerical 

analytic 

numerical 

analytic 

numerical 

analytic 

.0198 

1.394E5 

1.394E5 

4.334E-5 

4.382E-5 

6.942E-6 

7.044E-6 

.1147 

1.692E4 

1.692E4 

4.334E-5 

4 .381 E— 5 

6.942E-6 

7.044E-6 

1.090 

9.217E2 

9.217E2 

4.333E-5 

4.379E-5 

6.942E-6 

7.043E-6 

10.07 

1 .062E1 

1.062E1 

4.321E-5 

4.360E-5 

6.932E-6 

7 .027E-6 

100.1 

9.310E-3 

9.310E-3 

3.699E-5 

3.718E-5 

6.400E-6 

6.478E-6 

1059. 

5.089E-6 

5.089E-6 

2.014E-6 

2.019E-6 

8.741E-7 

8.799E-7 

10490. 

2.970E-8 

2.970E-8 

1.833E-8 

1.833E-8 

1.132E-8 

1.132E-8 
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TABLE 2 


Comparison of the Benchmark Numerical Simulation to the Analytic 
Solution for Secondary Manganese Ions as a Function of Ion Energy 
and Depth into the Aluminum Absorber. 


E ,MeV/amu 

$ Mn (10g/cm?E) 

*Mn (20g/cn1 » E) 

numerical 

analytic 

numeri cal 

analytic 

.0198 

1.772E-6 

1.780E-6 

5 . 704E-7 

5.768E-7 

.1147 

1 .772E-6 

1 . 780E-6 

5.704E-7 

5.768E-7 

1 .090 

1.772E-6 

1 . 7 7 9 E - 6 

5.704E-7 

5.767E-7 

10.07 

1 .767E-6 


5.696E-7 

5 .753E-7 

100.1 

1.504E-6 

1 . 503E-6 

5.242E-7 

5.291E-7 

1059. 

7 .797 E - 8 

7.806E-8 

6.880E-8 

6.918E-7 

10490. 

7.004E-10 

7.004E-10 

8.728E-10 

8.728E-10 













